Conductivity of quasiperiodic systems : a numerical study 
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We develop a new real-space method which allows to evaluate the Kubo- Greenwood formula for 
dc-conductivity of independent electrons in a static potential. We apply it to a numerical study 
of propagation modes in 3 dimensional quasiperiodic systems. These modes are strikingly different 
from those of periodic ones, with regard to the effect of disorder. In particular for Fermi energies in 
pseudo-gaps the conductivity can be stable or can even increase when disorder increases. 



Recently, studies on transport in general topologically 
disordered media have emphasized the importance to un- 
derstand how localisation of states and anomalous quan- 
tum transport are related |Q. For example, in disordered 
systems, close to a Metal-Insulator transition (at or close 
to mobility edges), localisation aspects of wavefunctions 
can be described by multifractal analysis [|| and related 
to critical exponents of conductivity || . Another exam- 
ple is the problem of localisation and transport properties 
in quasiperiodic systems. This has attracted a growing 
attention since the discovery of quasicrystals in 1984 j|. 
Historically, Kohmoto et al. (|] were the first to propose 
that, at least in ID systems, quasiperiodicity induced long 
range correlations giving rise to an intermediate state of 
localisation, namely "critical states" associated to sin- 
gular continuous spectra. Numerical studies of quantum 
dynamics of wave packets Q and Landauer resistances |t]] 
manifested peculiar features such as bounded resistance 
Pn < cste N a ( N,E >, for energies E in the spectrum of the 
Hamiltonian, with a(N, E) encoding the memory of these 
"quasiperiodic correlations" . More recently, correlations 
between quantum dynamics and localisation properties in 
Fibonacci chains have been rigorously analysed by means 
of renormalisation group treatments j^] and through new 
numerical methods (Iterated function systems @]) 

Experimentally quasicrystalline phases have unique 
electronic properties. They are characterized by a low 
conductivity, which increases when temperature or disor- 
der increases, and a proximity to a Metal-Insulator tran- 
sition JlOfl . Since real systems allways contain some de- 
fects either static, due to chemical or structural disorder, 
or dynamic due to phonons, it is of great interest to know 
how conductivity is affected by disorder in a quasiperi- 
odic system. 

In this context, bandstructure calculations of periodic 
approximants have been performed. They predict flat 



bands which are associated to very low Fermi velocity 
. These flat bands are associated to states that have a 
multifractal character. lljl^]. Estimates of the conduc- 
tivity within the Bloch-Boltzmann theory (and the relax- 
ation time approximation) lead also to small conductiv- 
ities. However the application of the Bloch-Boltzmann 
theory to these systems has been criticized ||lj|[ll). In- 
deed the propagation of electrons in the perfect quasiperi- 
odic structure is neither ballistic as in the case of periodic 
systems nor diffusive as in the case of disordered systems. 
One expects rather scaling laws of the form L(t) = At 13 
for the extension L(t) of a wavepacket, where (3 depends 
on the energy and on the Hamiltonian parameters. Thus 
in a first approximation the diffusivity will be given by 
D = L(t) 2 /3t = St 2 ' 3-1 where r is the finite lifetime 
induced by disorder. In the context of a model this has 
been confirmed by the mathematical work of Bellissard 
and coworkers [jl4j. Since states tend to be localised, it 
has also been proposed that the mechanism of conductiv- 
ity can be a hopping mechanism. In that case, according 
to the arguments Jl3|,[llJ , inelastic or even elastic scatter- 
ing could lead to an increase of the conductivity. 

In order to go beyond the Bloch-Boltzmann descrip- 
tion, and be able to test the various schemes that have 
been proposed, there is a natural starting point given 
by the linear response theory. This treatment does not 
make any assumption on the transport mechanism. In 
this Letter we present the first study, to our knowledge, of 
Kubo-Greenwood conductivity for independent electrons 
in static quasiperiodic potential at T — OA". Using a new 
real-space method, it will be shown that the electronic 
conduction in quasiperiodic systems differs strinkingly 
from the prediction of a Bloch-Boltzmann approach. In 
particular if Vdis measures the amplitude of the static dis- 
ordered potential we will show that the conductivity does 
not vary like a = o~o/Vdi S 2 (where <jq is independent of 
disorder) in the weak scattering limit. Furthermore the 
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variation of conductivity with disorder is rather complex, 
and in particular it depends strongly on the position of 
the Fermienergy with respect to pseudo-gaps. 

For independent electrons in a given static poten- 
tial at T = OK, the starting point of the method will be 
the following form of Kubo-Greenwood (E is the Fermi 
energy) : 



<tdc{E) 



2he 2 TT 

n 



Tr[V x 8{E-H) % S(E-H)] 



V x is component of the velocity operator along direc- 
tion x. 6(E — 7i) is the projector on eigenstates H of en- 
ergy E. Q, is the volume of the system. The factor 2 comes 
from the spin degeneracy. In this work, the evaluation 
of the Kubo-Greenwood formula is made through a new 
real space method. In a previous work Jl5[ we developped 
a method for calculating dc and ac conductivity using 
the formalism of orthogonal polynomials. Here, based on 
the same formalism, we developp a method which gives 
access only to dc-conductivity. In counterpart this al- 
lows an important gain in precision and numerical sta- 
bility. For convenience, after some simple algebra one 
can rewrite anc(E) as : 

fTDc(E) = lim F(E,t) 

t — >-oo 



F(E, t) = ^Tr [i (X(t) - X) 8{E - H) (X(t) - X)] 

Here X(t) = e lWt / s Xe~ lMf/s where X is the component 
along direction x of the position operator. At this stage, 
we define for each orbital j > : 



| $j(t) >=Xe- mt ' n | j > 

Here | $j (t) > is a normalised state. The calculation of 
conducivity can be reduced to : 



2he 2 n 



^V j (t)x<% j (t)\8(E-'H)\9 j {t)> 



For each initial state at | j >=| ^(4 = 0) >, we can 
take the origin at j. This means that < j \ X \ j >= 0. 
Then T>j (t) is given by 



<* J (t)\X 2 \t> (t)> 



The conductivity can thus be calculated from T>j(t) 
and the spectral quantity related to | 3>j(t) >. The 
spectral quantity is calculated by the recursion method 



]16[ . The original part of our method is the evaluation of 
| &j(t) >. Our method, to solve the Schrodinger equa- 
tion, avoids to use Runge-Kutta resolution algorithms, 
or numerical procedures of diagonalisation. It is based 
on the development of 5(E — Ti) already used in [ El| . 
which implies also a development of the unitary operator 
e -iHt/h Qn a ]-, as jg f orthogonal polynomials. Hereafter, 
we will take Chebyshev polynomials of the first kind as- 



sociated to the weigth p{E) = 1/(^^/46^ — (E — aoo) 2 ) 
and defined via the recursive relations : 



h(E)= 1 Qi(E) = 

'E-a, 



)n+i{E) = 



E — aoo 

2&oo 

n (E) - Qn-l(E) 



where and b^, are chosen as band parameters in ac- 
cordance with the ones of the true density of states (via 
a first recursion process). Then the vector under study 
reads : 



e n | j 



> =J2 hn (J dEp(E)Q n (E)e' 1 ^)Q n (n) \ j > 



ho = 1 and for other n h n — 1/2. As usual with or- 
thogonal polynomials the Q n (H) | j > are evaluated via 
the recurrence property [ jl6| . Amplitudes of | j > 

on this basis have a rather simple form connected to 
Bessel functions. The consequent interest is to get nice 
asymptotic behaviors for these coefficients which con- 
verge very quickly as n — > oo, given that lim„^oo J n {z) ~ 
-^== ( ) . It is this development that makes the cal- 
culation quick and precise. 

For each j the calculation is performed on a cube, cen- 
tered on the site j, which length is 100 sites (the cube 
contains about 10 6 sites). Also, for pratical calculations, 
the sum is performed over ~ 100 sites of origin. We find 
that this is sufficient for our purpose. We have performed 
numerous test (more details can be found elsewhere [fl7|) 
to ensure that the convergence is achieved. We estimate 
that the conductivity is calculated with an energy reso- 
lution of a few percent of the bandwidth which is enough 
for our purpose. As for the calculation of conductivity 
the total computing time is of the order of 100 hours on a 
HP735 for each conductivity curve cdc{E) (see below). 

We performed calculations of quantum diffusion and 
Kubo-Greenwood conductivity at T = OK. We consider 
a s-band tight-binding model on a simple cubic lattice 
with nearest-neighbor hopping which has already been 
studied by several authors (ll],|Hj. The hopping inte- 
gral is the energy unit (t = 1) and the on-site energies 



are given by £j 
random number 



e dls with e di , 



e [-^,+^] and e ja = ±V q 
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constraints to quasiperiodic correlations (Fibonacci se- 
quence) . This model allows a direct comparison between 
a quasiperiodic system and a periodic one since for V qp = 
one recovers the classical Anderson model with diagonal 
disorder. We note also that for Vdis = the Hamiltonian 
is separable. An eigenstate $(x, y, z) can be written as 
the product of eigenstates of the chains along each direc- 
tion &(x,y, z) — ^i(x) x ^2(2/) x ^3(2) the energy being 
the sum of the three energies E = E\ + E2 + E3. This 
means also that for a state that is initially localised on a 
site j one has ^j{x, y, z, t) — ^i{x, t) x ^(y, t) x ^ 3 (z, t) 
with ^i(x,t = 0) = 5(x — Xj) and similarly for y and 
z. Obviously the Hamiltonian is no more separable when 
Vdis is non zero. 

We study the quantum diffusion through Dj(t). For 
Vdis = the above relation shows that Dj (t) is the same 
as for a one dimensional model ||. In this study we fo- 
cus on the effect of disorder. Our results (see figure 1) 
, clearly show that conduction modes undergoes a tran- 
sition from non-ballistic to diffusive regime (i.e. Dj(t) 
is independent of t at large t). When the disorder in- 
creases the transition to the diffusive regime occurs at 
shorter times and the asymptotic value of Dj (t) tends to 
decrease. We note also that the fluctuations of Dj(t) are 
less important when the disorder increases. 

We studied also the conductivity <Jdc{E) and its vari- 
ation with the strength of disorder. In a metal one ex- 
pects a law of the form a = o^jVdis 2 (00 is independent 
of disorder). Indeed this is what we find for V qp = 0. 
For sufficiently small values of disorder our numerical re- 
sults are in good agreement with the prediction of the 
Bloch-Boltzmann approximation fl7| . However, the en- 
hancement of disorder in a quasiperiodic system leads 
to a different law. For instance in the figure 2 the con- 
ductivity decreases when the disorder increases but does 
not follow the law a = (Jo/Vdi S 2 (note that the density of 
states is nearly independent of disorder in this parameter 
range) . Instead we get approximately for each energy E a 
law of the form a = o-$/Vdi S a ' ■ Our results do not strictly 
follow a power law, probably because such laws apply 
only to the limit of infintely small disorder, and also due 
to the finite accuracy of our method. Depending on the 
energy E we find a = 0.4 — 0.8 in our calculation. 

In figure 3, we show the conductivity for a quasiperi- 
odic modulation V qp = 1.1. A striking result is that there 
are particular zones, identified by pseudo-gaps which 
seem quite insensitive to a tremendous increase of dis- 
order (a factor of 16 for Vdis 2 )- The conductivity varies 
monotonically for intermediate values of disorder not 
shown here. The inset in figure 3 shows the average den- 
sity of states. The density of states increases slightly with 
disorder in the pseudo-gaps. This increase compensates 
for the decrease of diffusivity leading to a nearly constant 
conductivity. In regions of high density of states, which 



correspond also to higher conductivity, both density and 
diffusivity decrease. 

In figure 4 we show aoc(E) for V qp — 2.5 and different 
values of Vdis (again the conductivity varies monotoni- 
cally for intermediate values of disorder not shown here). 
The variation of conductivity with energy is rather com- 
plex. We note however that the conductivity increases 
with Vdis for some energies that correspond to pseudo- 
gaps. In contrast to the previous cases there is an impor- 
tant change of the density of states. It increases where 
the conductivity increases. It means that the electronic 
structure is deeply modified by the disorder. Thus we 
prove that the localization induced by a quasiperiodic 
potential can be destroyed by a disordered potential. 

A direct comparison with experimental results on qua- 
sicrystals is difficult, since our model does not provide a 
realistic description of the electronic structure and of the 
local atomic order. However, quite remarquably, when 
the Fermi energy lies in a pseudo-gap, which is the case 
experimentally, the behavior of the conductivity is remi- 
niscent of the experimental observation that the conduc- 
tivity increases with disorder or with temperature ]l0| . 

To conclude, the contribution of this paper is 
twofold. We have developped a new real-space method 
which allows to study electronic propagation and to cal- 
culate the conductivity of quasiperiodic systems (through 
the Kubo- Greenwood formula). This opens new ways of 
investigating the transport properties of these systems, 
but of course the method is not restricted to quasiperi- 
odic potentials. Furthermore our first numerical study of 
conductivity demonstrates the complexity of electronic 
transport in quasiperiodic systems with disorder. In 
the context of anomalous localisation, it shows how a 
quasiperiodic potential can affect quantum transport in 
a specific way, when compared to a periodic or a dis- 
ordered potential. For sufficiently small disorder and 
quasiperiodic potential the reduction of conductivity is in 
qualitative agreement with the prediction of scaling laws. 
However a particularly interesting result is the inhomo- 
geneous variation of conductivity with Fermi energy. In 
zones of pseudo-gaps the conductivity can be stable upon 
enhancement of disorder or can even increase with disor- 
der. 

ACKNOWLEDGEMENTS 

This work was partially supported by a NATO grant 
CRG 941028 which is gratefully acknowledged. 

FIGURE CAPTIONS 

Figure 1 : T>j (t) represented for several initial sites j 
and V qp = 0.9. Time is in units of 2h/W where W is the 
bandwidth. Lengthes are in unit of the nearest-neighbors 
distance. (a,)Vdi S = 2 (b) Vdis — 2y/2. The thick line is 
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an aid to visualize one of the curves. 

Figure 2 : Conductivity a D c(E) for V qp = 0.7 (av- 
erage value over sites "j" of lim^oo Dj (t) Nj (E) where 
Nj(E) is the spectral weight on the state | <&j(t) >). The 
energy unit is the hopping integral t. Inset : average 
density of states, (a) V dis = 2. (b) V dis = 2\/2 

Figure 3 : Conductivity a D c{E) for V qp = 1.1 as a 
function of the Fermi energy E for different values of the 
disorder parameter. Inset : average density of states, (a) 
V dis = 1/V2 (b) V dis = 2V2. 

Figure 4 : Conductivity rr DC {E) for V qp = 2.5 as a 
function of the Fermi energy E for different values of the 
disorder parameter. Inset : average density of states. 
(a,)V dis = 2 (b) V dts = 4V2. 
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